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pression for the hypermultiplets' leading-order contribution to the free energy, and in the 
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hypermultiplets' leading-order contribution to the free energy at zero mass. 
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1. Introduction and Review 

Understanding the phase diagram of Quantum Chromodynamics (QCD) is essential for un- 
derstanding a variety of physical systems, such as the early universe and neutron stars. Much 
of the interesting thermodynamics of QCD, such as the deconfinement transition, occurs 
when traditional perturbative methods are unreliable. The lattice formulation of the the- 
ory circumvents this problem by direct calculation of the path integral. Implementing a 
baryon number chemical potential on the lattice is difficult, however. We are thus led to 
study the thermodynamics of model theories similar to QCD, but that offer more tractable 
strong-coupling analysis. In particular, the AdS/CFT correspondence, and more generally 
gauge-gravity duality, permits quantitative study of the thermodynamics of strongly-coupled 
non-Abelian gauge theories that share some features with QCD. 

In this paper we will use the AdS/CFT correspondence to study the = 4 supersym- 
metric SU{Nc) Yang-Mills (SYM) theory, in the 't Hooft limit of large-A^c, when the 't Hooft 
coupling is large, A = gyj^Nc ^ 1. The M = 4 SYM theory is very different from QCD at 
zero temperature, being a conformal field theory (CFT) and hence lacking confinement, a 
mass gap, and chiral symmetry breaking. At finite temperature, however, the AA = 4 SYM 
theory resembles the high-temperature phase of QCD, in particular by exhibiting screening 
of color charges. The AA = 4 SYM theory at finite temperature has no phase transitions, 
however. Being a CFT, the theory has no intrinsic scale to set a transition temperature, so 
no transitions can occur. 
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The AdS/CFT correspondence is the conjecture that the AA = 4 SYM theory, in the lim- 
its described above, is equivalent to type IIB supergravity formulated on the background 
spacetime AdS^ x S^, where AdS^ is five-dimensional anti-de Sitter space and is a five- 
sphere [1—3]. This spacetime arises as the near- horizon geometry of Nc ^ oo coincident 
D3-branes in type IIB string theory. Thermal equilibrium of the A/" = 4 SYM theory cor- 
responds to non-extremal D3-branes, whose near-horizon geometry is five-dimensional AdS- 
Schwarzschild times [4, 5]. The temperature T of the A/" = 4 SYM theory is identified with 
the Hawking temperature of the AdS-Schwarzschild black hole. 

The A/" = 4 SYM theory, unlike QCD, only describes fields in the adjoint representation 
of the gauge group. We will introduce fields in the fundamental representation in the form 
of A^=2 supersymmetric hypermultiplet "flavor" flclds. We will introduce a finite number 
Nf of them, so that Nf <C Nc in the 't Hooft limit, and work to leading order in Nf/N^. 
To this order, the beta function of the theory is zero and hence the theory remains (approx- 
imately) conformal. We may introduce an A/" = 2 supersymmetry-preserving mass m for the 
hypermultiplets, however, which of course explicitly breaks the conformal invariance. 

The Nf hypermultiplets possess a U{Nf) global symmetry analogous to the vector flavor 
symmetry of QCD. As in QCD, we identify the U{1)b subgroup of U{Nf) as baryon number. 
We will introduce a finite U{1)b density (J*), where is the U{1)b current, or equivalently 
we work with a finite chemical potential jj,. 

The flavor fields appear in the supergravity description as Nf D7-branes [6]. In the 
Nf ^ Nc limit, the contribution that these D7-branes make to the stress-energy tensor is 
dwarfed by the contribution of the Nc D3-brancs. To leading order in Nf /Nc, then, we may 
neglect the back-reaction of the D7-branes on the geometry: they are probes. The D7-brane 
action is then the Dirac-Born-Infeld (DBI) action. The hypermultiplet mass m appears in 
the supergravity description as the asymptotic separation of the D3-branes and D7-branes 
in a mutually orthogonal direction. The U{Nf) flavor symmetry appears in the supergravity 
description as the U{Nf) worldvolume gauge invariance of the Nf D7-branes. The U{1)b 
current is dual to the U{1) worldvolume gauge field A^. To obtain a finite baryon number 
density (J*) in the SYM theory, we nuist introduce a nonzero gauge field time component, 
At{r), where r is the radial coordinate of AdS. 

The SYM theory has three dimensionful parameters: the temperature, T, the hypermul- 
tiplet mass m, and the baryon number density, (J*). We will use the D7-brane description 
of the flavor fields to derive exact expressions for their free energy in two limits: T = and 
m = 0. To place our results in context, we will now review the known thermodynamics of 
the fiavor fields in the SYM theory, as deduced from supergravity calculations. 

Consider first zero temperature, zero mass, and zero density in the SYM theory. The 
D3-branes and D7-branes then have zero asymptotic separation. The D7-branes are extended 
along AdS^ x inside of AdS5 x S^, where the worldvolume is the maximum- volume 
equatorial C [6]. We first turn on m, which means finite asymptotic separation of 
the D3-branes and D7-branes. The D7-brane then has a non-trivial embedding in the near- 
horizon geometry. The D7-brane wraps the equatorial C at the boundary of AdS^. 
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Moving into the bulk of AdS^, the position of the can change, that is, the D7-brane "shps" 
on the S^. The volume of the shrinks as we move into AdS^ and indeed may ultimately 
vanish at some radial position r', as allowed by topology (the is a trivial cycle). The 
D7-brane then does not extend past r': it appears to "end" at r' [6]. The worldvolume scalar 
that describes the embedding is dual holographically to an operator Om given by taking ^ 
of the SYM Lagrangian. In particular, Om includes the mass operator of the fermions in the 
hypermultiplet as well as m times the mass operator of the scalars, as well as couplings to 
adjoint scalars. The exact operator is written in ref. [7]. For our purposes, thinking of Om as 
the mass operator will be sufficient. 

Now turn on the temperature T. The spacetime is now AdS-Schwarzschild, and a horizon 
is present at some radial position ru- Two classes of embedding are now possible. The first 
are the finite-temperature analogues of the embeddings described above: D7-branes that end 
outside the horizon, i.e. r' > rn- These are called "Minkowski" embeddings. The second 
class are D7-branes in which the C shrinks as we move away from the boundary, 
but never collapses to zero volume. The D7-brane then extends to, and intersects, the AdS- 
Schwarzschild horizon. These are called "black hole" embeddings. A topology-changing phase 
transition exists between these two types of embeddings, which appears in the SYM theory 
as a first-order phase transition associated with the flavor fields as we change m/T [8-15]. 

Now turn on the baryon number density (J*)- The U{1) gauge field on the D7-brane 
worldvolume now has nonzero At{r). The asymptotic value ^t(oo) gives the value of /x. In 
the canonical ensemble, the topology-changing transition observed at zero density persists for 
small density and then ends in a critical point [7]. A region of instability exists, with the line 
of first-order transitions boundary [7]. 

A claim was made in ref. [7] that when Ftr is nonzero only black hole embeddings are 
physically allowed. Minkowski embeddings were dismissed as unphysical because, if the D7- 
brane ends outside the horizon, the radial electric field lines Ftr have no place to end. If a 
source, a density (J*) of fundamental strings stretching from the D7-brane to the horizon, 
is introduced to accommodate the field lines, the tension of the strings will overcome the 
tension of the D7-brane, drawing the D7-brane into the horizon and producing a black hole 
embedding. Minkowski embeddings were thus excluded from the analysis of ref. [7]. 

In refs. [16,17], however, an argument was made that Minkowski embeddings were phys- 
ical and should be included. One of the main objections to the claim of ref. [7] was that the 
black hole configurations cannot realize chemical potentials below a critical value, which at 
zero temperature wc will show is fi = m. What we find is that this "incompleteness problem" 
is completely resolved by including Minkowski branes with a constant At{r), that is, with 
non- vanishing chemical potential but zero density^. The authors of refs. [16,17] argue that 
At is not a gauge-invariant parameter, but in the thermodynamic setting where spacetime 
has a compact Euclidean time direction, the Wilson line of At is gauge invariant and hence 
physical. The constant part in At{r) is the chemical potential according to the standard rules 

^While this paper was in preparation ref. [18] appeared, wliere the same conclusion was reached. 
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of AdS/CFT, since it represents the leading near-boundary behavior. For any configuration 
that satisfies the boundary condition that Aj(r) vanish in the infrared (deep in AdS), this 
definition becomes equivalent to the description fJ- = J F^t advertised in refs. [16, 17]. The 
latter definition misses the important case that At{r) is simply constant, however, which is 
always a solution in the absence of sources. 

We will find analytic solutions for the zero-temperature limit of black hole embeddings, 
and will argue that these are thermodynamically favored when fi > m and (J*) > 0. This is in 
complete agreement with the claims of ref. [7], where only (J*) > was considered. When fj, < 
m, Minkowski embeddings with constant At{r), and hence (J*) = 0, are thermodynamically 
favored. The transition between the two is second order. Indeed, when /i > m we see, as 
expected, a free energy whose form is consistent with spontaneous symmetry breaking for a 
scalar field with a wine-bottle potential. When ^ = m the minimum of the scalar's potential 
returns to the origin. Note in particular that we will not see a Fermi surface: the system 
prefers to condense scalars. 

We can also generalize our zero-temperature analysis to include nonzero currents (J*), for 
spatial index i. At zero temperature, these currents will not dissipate and can be introduced 
as a constant background. We can generalize further, from a D7-brane probe to a D5-brane 
probe [19,20]. Indeed, our analysis should generalize to many probe D-brane systems. 

At finite temperature, we are also able to provide an analytic expression for the hyper- 
multiplets' contribution to the free energy at leading order in Nj/Nc in the limit m = 0. The 
result exhibits no interesting phase structure, in accord with the phase diagram of ref. [7]. 

Another holographic model of QCD at finite baryon density, the D4/D8 system, and 
specifically the Sakai-Sugimoto model [21], has received a great deal of attention recently 
[22-29]. In this case, the supergravity background is provided by D4-branes and the flavor 
branes are D8-branes and anti-D8-branes. As in ref. [30], in the course of our analysis we will 
flnd a D7-brane configuration similar to the D8-brane in the Sakai-Sugimoto model, but, also 
as in ref. [30] , we will find that these configurations are always thermodynamically disfavored. 
Indeed, these solutions should not be compared to any of our other solutions, as they are non- 
supersymmetric even at zero density and describe a SYM Lagrangian distinctly different from 
that of the SYM theory described above. 

Most of the work on the Sakai-Sugimoto model at finite density is based on analogues 
of Minkowski embeddings, with strings attached (or alternatively a wrapped D4-brane pro- 
viding the required source). Our results, as well as the earlier work of ref. [7], show that 
in the analytically tractable D3/D7 system, such configurations are not solutions of the full 
source-plus-brane system, and do not even provide a good qualitative picture of the phase 
structure. This at least suggests that any analysis based on approximating a full solution to 
the non-abelian DBI action relevant in the Sakai-Sugimoto model with a source-plus-brane 
configuration has to be approached with some caution. Presumably a solution like the smooth 
black hole embedding for the D7-brane discussed here exists in that case as well, but it most 
likely will involve non-abelian configurations of scalars and gauge fields as well as the tachyon. 

This paper is organized as follows. In section ^ we present our zero-temperature analysis. 
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In section ^ we present our zero mass analysis. We conclude in section We collect various 
necessary integrals in the Appendix. 

2. The Limit of Zero Temperature 
2.1 Analytic Brane Embeddings 

In type IIB string theory, we consider a system of Nc D3-branes and Nf D7-branes aligned 
in flat ten-dimensional space as 

Xq Xi X2 X3 X4 ^5 Xq Xj Xg Xg 

D3 X X X X (2.1) 

D7 xxxxxxxx 

If we separate the D3-branes and D7-branes in a mutually orthogonal direction, in the rcg- 
xg plane, an open string stretched between them will have a mass equal to the separation 
times the string tension. This mass appears as the hypermultiplet mass m in the D3-brane 
world volume theory. We take the usual AdS/CFT limit, Nc ^ 00, gs ^ with QgNc fixed 
and gsNc ^1 [1]. We obtain the near-horizon geometry of D3-branes, AdS^ x S^, with the 
Nf < Nc D7-branes extended along AdS^ x S^. 

The DBI action describing the embedding of any probe D-brane is non-linear. Even when 
the symmetries are strong enough to restrict all the embedding functions to depend on a single 
variable only, the resulting ordinary differential equation of motion can often only be solved 
numerically. Sometimes symmetries provide integrals of motion that permit an exact solution 
of the system. For example, in the Sakai-Sugimoto model [21] translation invariance along 
the direction x transverse to the D8-branes ensures that the D8-brane action only depends 
on derivatives of x and not on x itself. We will show that the same is true for the D3/D7 
system at zero temperature even in the presence of certain worldvolume gauge fields. From 
the SYM perspective, this is the statement that we may determine the free energy exactly 
for non- vanishing currents J^, and in particular for finite density (J*). 

In the original analysis of ref. [6], the embedding of the D7-brane in AdSs x was 
characterized by a slipping mode 6{rQ). That is, the AdSs x metric was written as 

ds^ = -l^T^ dx^'dx" + ^drj + (de^ + sin^ 6 dsli + cos^ 6 dsls ) (2.2) 

where = 0,1,2,3 denote the four worldvolume directions of the D3-brane, ry^jy is the 
(3-|-l)-dimensional Minkowski metric, R is the curvature radius^, ds^i and ds^a are the 
metrics of unit 1- and 3-spheres, respectively, and rg is the radial direction in the full space 
transverse to the D3-branes. In terms of the ten-dimensional space in which the A^c-D3 branes 
are embedded, r| = X^a-^a' where a runs from 4 to 9. The action of a D7-brane wrapping 

^Starting now, we work in units with R = 1. In particular, we convert between supergravity and SYM 
quantities with a' = X~^^^. 
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AdSs X S"^ is proportional to the volume of the S^, which goes as cos^0, so we will find a 
non-trivial potential for the embedding scalar ^(rg) in this case, and no simple integral of 
motion is apparent. However, if we write the background metric as 



■'^dsls + dy^ + dz^) 



(2.3) 



a significant simplification occurs. In this case, y{r) and z{r), which stand for xg and xg, 
are the scalars describing the embedding, and without loss of generality we can use the U{1) 
symmetry that rotates them into one another to set z = 0. r is a radial coordinate in four of 
the six transverse directions, so we have rg = + + z^. All the y-dependence of the action 
appears via rg in the warp factors ^(rg). For the embeddings we are interested in, with a D7- 
brane wrapping the four directions and four of the six transverse directions, all the warp 
factors will drop out from the action, and we are effectively solving for D7-brane embeddings 
in flat space. The same simplification occurs in all supersymmetric probe brane embeddings 
where background and probe branes have four relative ND directions. The action now only 
depends on derivatives of y, and the resulting constant of motion is sufficient to integrate 
analytically the equations of motion. This cancellation of warp factors is spoiled at finite 
temperature or in the presence of electric or magnetic background fields in the SYM theory. 
This coordinate system was originally used in ref. [31] to find the exact meson spectrum of 
the SYM theory from fluctuations of the D7-brane. More recently, it was used in ref. [30] 
to flnd new non-supersymmetric D7/anti-D7-brane embeddings. To translate back to the 
original 9{r) embedding scalar, we can always combine the deflnitions y = rgsinS and (at 
z = 0) Tg = + to obtain tan 9 = 
The action for the D7-branes is 



where Toy is the D7-brane tension, are worldvolume coordinates, gab is the induced world- 
volume metric and Fab is the worldvolume U{1) gauge fleld. Introducing the fleld strength^ 
Frt = A[{r) and the embedding scalar y{r), the zero-temperature D7-brane Lagrangian is 




(2.4) 




(2.5) 



where M = NjTdjVs with V3 = 27r^ the volume of a unit S^. In terms of SYM theory 
quantities, 




We will divide both sides of eq. by the volume of M^'^ and henceforth work with action 

(and in a moment, free energy) densities, so that Sd7 = / drL. 



Si 



Throughout this paper, we work in a gauge with Ar — 0. 
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Thanks to the warp factor cancellation, only derivatives of y(r) and At{r) appear in the 
Lagrangian, and we find two conserved charges 



^ = y' ^ -c ^^^3 ^ ^ 

(2.7) 

The ratio implies 

We can algebraically solve for y'{r) and A[{r) in terms of the integration constants c and d, 

V = , A', = "/'^'-"''^ (2,9) 

These can be integrated using incomplete Beta functions (whose properties we review in the 
Appendix). The result depends on the sign of (^2va')'^ ~ When c = d = 0, we obtain the 
solution with y'{r) = and A[{r) = 0, so y{r) and At{r) are constants. If At{r) = 0, we 
recover the original embedding of ref. [6] for a single straight D7-brane. This embedding is 
depicted in figure (Q a.). Solutions with ^2ncx')'^ ~ positive and negative are also depicted 
in figure (||). We address each of these in detail in the following sections. 
The action evaluated on these solutions is 



fA 

Sd7 = -M I drr^ 



(2.10) 



■'^ ' ^ (2™')^ 



The lower endpoint of integration will depend on the sign of (27ra')^ —(?■ The integral diverges 
if we integrate to r = oo, so we regulate the integral with a cutoff at r = A. The divergence 
is present for d = c = 0, i.e. for straight D7-branes with At{r) = 0, 



^A -1 

Sq = -M drr^ = --Mh!^ (2.11) 
JO 4 



In what follows we will obtain the renormalized on-shell action, Sren, by subtracting this 
divergence, 

Sren = Hm {Sm - Sq) (2.12) 

A— »oo 

In the SYM theory, the thermodynamic potential of the grand canonical ensemble 
is given by O = —Sren- We Legendre transform to obtain the free energy density in the 
canonical ensemble, F = O + jid. Recall that at zero temperature, free energy and energy are 
the same, so at zero temperature F is also the energy density. The conserved charges c and 
d determine {Om) and (J*) as follows. 
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(a.) (b.) (c.) 

Figure 1: Cartoons of the allowed D7-brane embeddings at zero temperature. The vertical axis is 
one of the directions transverse to the D3-branes but parallel to the D7-branes. The horizontal axis 
is an overall transverse direction, (a.) A Minkowski embedding corresponding to y(r) being constant. 
(^■) — < solution with turnaround describing a brane/anti-brane embedding as in 

ref. [30]. (c.) The zero-temperature — > horizon-crossing "black hole" embedding. 



where in each case when we vary one field we hold the other fixed. We then have 

6Sd7 = j dr (^-^^^dr5At{r) + JI±-^dr6y{r)^ = d5At{^) - c6y{^) (2.14) 

where we demand that 5At[r) and 5y{r) are always zero at the lower endpoint of the r 
integration. If we vary At{r) while holding y[r) fixed {6y{r) = 0), we find (J*) = d, and 
similarly we find {Om) = (27ra')c. Starting now, we will stick to the notation c and d and 
refer to these as the condensate and density. 

Using f2 = —Sreni we can determine whether solutions with nonzero d and c are thermo- 
dynamically favored relative to the single straight D7-brane with zero gauge field. Sren > 
means the configuration is favored and Sren < means the configuration is disfavored. 

2.2 Brane/Anti-Brane Embeddings 

We begin with ^2na'y^ — < 0. In particular, we will begin with d = solutions, as analyzed 
recently in ref. [30]. For such d = embeddings, y'{r) diverges as at a critical radius 

ro = M~^^^c^^^ . This behavior indicates that the D7-brane has a turn-around point where the 
D7-brane smoothly matches onto a second branch. This solution thus describes a D7/anti- 
D7-brane pair connected by a smooth neck, as depicted schematically in figure (|| b.). 
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In the asymptotic region, we may specify the distances of the D7-brane and the anti- 
D7-brane to the D3-branes separately. In other words, we may associate two masses to such 
a configuration. We denote these as m and m. A D7-brane solution exists for any m and 
{Om)- The mass (or asymptotic separation) of the anti-D7-brane is then fixed in terms of 
these two input parameters. From tlic field theory point of view, we would like to interpret 
the two masses m and fh as input parameters, so that the above construction dynamically 
determines the value of the condensate (Cm)- Prom the supergravity picture, we can see that 
the condensate only depends on the difference of the masses, m — fh. This is another reflection 
of the fact that the embedding equation lost all knowledge of the warp factors, and effectively 
we are solving for flat-space embeddings. Only the relative position of the D7-branes and 
anti-D7-branes matters, not their distance from the origin. As noted in rcf. [30], however, all 
the connected D7/anti-D7-brane solutions have higher energy than the straight D7/anti-D7- 
brane pair, and hence will be meta-stable, at best. The true ground state of the system has 
vanishing {Om)- 

Now we introduce nonzero d. To connect the two halves of the D7-brane smoothly we 
need the gauge field and its derivative to be continuous at the neck. In particular, the field 
strength on the two branches should be equal in magnitude but opposite in direction. The 
baryon number density is then equal and opposite at the two ends, so this configuration does 
not describe any net baryon number. The Lagrangian has a global U{1) x U{1) symmetry, 
and, at least in the grand canonical ensemble, we can imagine turning on different chemical 
potentials for the two i7(l)'s at the level of the Lagrangian, even though the symmetry is 
spontaneously broken to the diagonal U{1)b in the state represented by the connected D7- 
brane configuration. If we wanted to have a branc/anti-brane system with finite baryon 
density wc would have to include an explicit source at the neck. Such additional sources will 
back-react and alter the configuration, so we will not discuss this option further. 

In any case, these brane/anti-brane systems correspond to a different Lagrangian than 
that for systems with a single D7-brane reaching the asymptotic region. The dual field 
theory for the brane/anti-brane is A/" = 4 SYM coupled to two hypermultiplets that preserve 
opposite J\f = 2 supersymmetries, so that the full theory is non-supersymmetric. When 
comparing all possible bulk configurations that correspond to hypermultiplets preserving A/" = 
2 super symmetry, the brane/anti-brane configurations never contribute. We will therefore 
only compare these solutions to the straight D7-brane and anti-D7-brane. 

For nonzero d, keeping in mind that for the brane/anti-brane system this does not cor- 
respond to a net baryon density, the D7-brane reaches the turn-around point 



The integrals for y(r), At{r) and the regulated on-shell action are performed in the Appendix. 




(2.15) 



We find 
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and At{r) is simply (27rl')^ c times y(r), as shown in eq. (^). Notice y{ro) = At{ro) = 
0. The mass and chemical potential are given by the asymptotic values at r — > oo, with 
lim,.^oo B |, = 0. The renormalized on-shell action for the D7-brane is 

S„„.-lA-/=(c^-p^)^%(-Hi) (2.17) 

The action for the anti-D7-brane is identical, so the total action is twice Sren- Notice that 
Sren < 0, SO these connected brane/anti-brane configurations always have a higher free energy 
than the straight brane and anti-brane with constant At{r), and hence are thermodynamically 
disfavored, just as in the d = case. 

At this point we may also consider c = / so that (2™^"*^^ ~ ^' ^^^i^ configuration 
shares a similar feature with the brane/anti-brane configurations: it does not satisfy the 

d2 



required UV boundary conditions for a theory with a single D7-brane added. When 



(2™') 



= 0, we have y'{r) ~ so the D7-brane misses the D3-branes at the origin and returns 
to the asymptotic region at r — > 00. We will therefore not study the c = case further. 

2.3 Finite-Density Black Hole Embeddings 

We now consider — > 0. In this section we will show that these are the zero- 

temperature limit of black hole embeddings. First, notice that these solutions have no turn- 
around point. They describe D7-branes alone. 

For a (^2na')'-^ — > solution, suppose we take c = 0. The solution with constant y{r) 
is still allowed. Geometrically, this is the zero-temperature Minkowski embedding, but now 
with nonzero A[{r). As mentioned before, the constant charge density along the D7-brane 
requires a source at r = 0, so this embedding is not physical unless we add extra strings 
connecting the D3-branes and D7-branes, or some other source. We will return to these c = 
embeddings in the next subsection. 

We can avoid this issue, however, by demanding that the D7-brane touch the horizon at 
r = 0. Since the horizon is located at = rg = -|- y^, this immediately implies a boundary 
condition on y 

y{r = 0) = 0. (2.18) 

The embedding satisfying this boundary condition is the zero temperature analogue of black 
hole embeddings. We will stick to this name, even though at zero temperature "horizon 
crossing" would be more appropriate. Using the change of variables tan^ = y/r, we see that 
the slope of y{r) at r = sets the value of 9 at the horizon. This black hole embedding is 
displayed schematically in figure (|| c). 

The integral for y'{r) in eq. (|2.9| ) is done in the Appendix, with the result 
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^ (27ra')^ 

and again, At{r) is j2m^c ti™6s y(r). Notice that y(0) = A((0) = 0. Identifying m and /x 
from the asymptotic values y{oo) and ^^(oo), where as r ^ oo the incomplete Beta function 
becomes B (g, |), we find 

c = 7AA(27ra )^ (/i^ - m^) m (2.20a) 
= 7AA(27ra')3 (/x^ - m^) /i (2.20b) 

ZTTCt 

where we have defined the constant 

"'^{l'' {hi)) p-^i* 

Note that since these black hole embeddings were valid only when — > 0, we 

can immediately see that they only exist for /i > m. In a theory with a mass gap, chemical 
potentials of magnitude less than the mass gap are basically trivial: they do not lead to any 
finite density. In the canonical ensemble they never appear, since zero density corresponds 
to fi = m. In the grand canonical ensemble we should of course be able to dial /i all the 
way down to zero, but the physics for < m is still expected to be trivial in this range, 
with d = and the free energy being independent of /i. These expectations will be verified 
when we study Minkowski embeddings in the next subsection. In this sense, the black hole 
embeddings cover the whole range of physically interesting values of //. 

The functional form of the condensate c(;U, m) is exactly that of a charged boson in the 
presence of a chemical potential. The chemical potential leads to a negative mass squared for 
the boson. Once this is larger than the positive mass squared in the Lagrangian, and if the 
scalar has a quartic coupling, the scalar will condense, its expectation value settling to the 
new minimum. Eqs. ( p. 20 ) can thus be viewed as determining the effective quartic potential 



at strong coupling. Notice that this also confirms the field theory expectation that at finite 
density the ground state is described by scalar condensates, not by a Fermi surface. 
The renormalized on-shell action in this case is 

Sren = l^-'/^M-^/' - 0^^^ = 1,^(2-0.')^ " m^f (2-22) 

Notice that Sren > 0, so these embeddings are thermodynamically favored relative to the 
straight D7-brane with no gauge field. In the next subsection we will compare to the straight 
D7-brane with a nonzero gauge field. Notice that we observe no instability in Vt. Stability 
requires ^ > 0, which is clearly satisfied in eq. ( |2.2C| ). This is consistent with the numerical 
results of ref. [7] in the zero-temperature limit. 

The thermodynamic potential in the canonical ensemble is 



SO that i*" > 0. We will see in the next subsection that the black hole embedding still turns 
out to be the configuration with lowest F, however. 

2.4 Finite-Density Minkowski Embeddings 

Let us return to c = and constant y{r). These are solutions everywhere away from r = 0. 
These are straight D7-branes, which do not obey the boundary condition y{0) = 0. In other 
words, when r = the D7-branes are still a distance y away from the D3-branes. The mass 
m is given by m = 

The integral for At{r) will be unchanged, however. In particular, we still integrate from 



to r = oo, and At{0) = 0. From eq. (2.20) we read off the relation between ^ and d 



d = 7AA(27ra')V^ (2-24) 
The integral for the on-shell action is also unchanged, so we simply set c = in the first 



equality of eq. (|2.22[) , 



Sren = M-^l\2Tia'r^l^ d^'"" = ^-^M {2^a'f (2.25) 

and hence = —Sren ~ ~/^^- Naively this for any nonzero d (or /x > m), is smaller (more 
negative) than the one for the black hole embedding, eq. ( ^.22] ). Notice, however, that unlike 
the nice interpretation we found for the relations among c, d, /i and m for the black hole 
embeddings, this result for $7, a pure quartic in fi, is completely independent of the mass and 
therefore appears to be unphysical. For example, in the decoupling limit m — > oo, O should 
go to zero. 

The problem of course is that the Minkowski embedding, as it stands, is not consistent: 
we completely ignored the contribution to O due to whatever object sources the D7-brane 
gauge field. The simplest source is a finite density d of fundamental strings stretching from the 
D3-branes to the D7-brane. These D7-branes with strings attached are not solutions to the 
full equations of motion following from the combined DBI and Nambu-Goto actions, however. 
Naively, we might think that the string back-reaction on the D7-brane can be neglected in 
the large- A'c and large-A limit. For a single string this is certainly true. The prefactor of the 
Nambu-Goto action scales as \/A, while from eq. (|2.6| ) we see that the prefactor M of the 
DBI action scales as NjNcX. We want a finite density d of strings, however. If we keep all 
geometric distances of order one in units of the AdS curvature radius, the mass will scale as 

so we want the chemical potential also to scale as \/A- From eq. ( |2.20[ ), we can see that 
the density d will then scale as NfNc\^- With this, d times the Nambu-Goto action has the 
same NfNcX scaling as the DBI action, and so the back-reaction is order one. 

To include the back-reaction one has to re-solve the equations of motion including the 
extra source term. We already found the most general solution to the equations of motion 
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with sources localized at r = 0, so the back-reacted solution including any such sources 
must be within this class. The only well-behaved solution in this class which asymptotically 
becomes a single brane is the black hole embedding, so it must be the back-reacted solution. 

To show that this is physically reasonable we want to show that the brane-plus-strings 
configuration, which could serve as a consistent initial data for a full time-dependent physical 
solution, has higher energy than the black hole embedding. We add a term dm to the action 
representing the finite density d of strings with length y = {2'Ka')m. The full free energy 
including the D7-branes and strings in this case is 

^^L^ = _/ + 4^3^ = /,3(4^_^) (2.26) 

which is positive at = m and is disfavored relative to the black hole embedding at all 
values of n > m. We can also allow some mixture, where the D7-brane embedding satisfies 
the boundary condition y{r = 0) = yo foi" some m > > 0, with a density d of strings 
extending from y = to y = yo. In figure (||), we imagine sliding the black hole embedding 
to the right a distance yo, such that the total asymptotic separation remains y, and then 
introducing the strings. The free energy is then 



, =-(f^'-(m-^Y)\ 4^ (^.' -(m-^Y)-^ (2.27) 

and we can easily check that in the relevant regime ^jl > m > > the derivative of this 
expression with respect to yo is strictly positive. The D7-brane with strings attached can 
continuously lower its energy until it turns into the black hole embedding with all strings 
dissolved in the D7-brane. 

While none of the finite-density Minkowski embeddings solve the equations of motion, 
the trivial Minkwoski embeddings with constant At{r) do. As only Frt enters the action 
and not At{r) itself, these embeddings have F = 0, = d = 0, and ;U is a free parameter. 
For n < m these are the only allowed configurations, so they dominate the ensemble. The 
fact that they are indistinguishable from the vacuum is completely natural from the field 
theory point of view: the lightest charge carriers have mass m, and for // < m no non-zero 
density can be produced. When fj, > m these solutions still exist, but in that regime the 
black hole embedding, with negative fi, dominates the ensemble. At /i = m we hence have a 
second order phase transition between the trivial Minkowski embedding and the black hole 
embeddings. Both Q and its first derivative with respect to fi are zero at that point, but the 
second derivative jumps from zero for the trivial Minkowski embedding to —2^J\f{27ra')^m? 
for the black hole embedding. Again, this is exactly what we expect for a bosonic condensate 
with a non-trivial effective quartic interaction. The critical exponents are given by their 
mean-field values, that is, for fi = m + e we have ~ e^, c ~ e, d ~ e. 

In the canonical ensemble, the trivial Minkowski embeddings do not even contribute 
because they have d = 0. The finite-density Minkowski embeddings with strings attached that 
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we considered above are not solutions of the equations of motion, so they cannot contribute. 
We thus conclude that, in the canonical ensemble, only black hole embeddings contribute. 
This analytically confirms the results of ref. [7], at least at zero temperature. 

Lastly, notice that the configuration we mentioned above of a single straight string stretch- 
ing between a D3-brane and a D7-brane in flat space is BPS, preserving half of the super- 
symmetry. We should ask why a density of such strings is not also BPS, and hence the 
dominant configuration? From a field theory point of view, the single string represents a 
single fundamental-representation "quark." As the force between colored charges is Coulom- 
bic, we may introduce a single quark into the system and still make the total state of the 
system gauge invariant by placing an anti-fundamental color charge at infinity. A finite den- 
sity of such states, as would be represented by the configuration we studied above, is not a 
physical configuration, however, because it has a finite color density. It does not represent 
a gauge invariant state that contributes to the partition function. Both at strong and weak 
coupling a scalar condensate forms at finite density, breaking the SU{Nc) gauge group to 
SU{Nc — 1) X U{1). In other words, the system moves out onto the Higgs branch. No BPS 
states exist on the Higgs branch. For example, a finite density of baryons is gauge invariant 
but is not BPS. Our analysis is consistent with this field theory picture: the D7-brane-with- 
strings configuration is not a supergravity solution, and in the SYM theory we indeed find 
that a scalar condensate forms. 

2.5 General Zero-temperature Probe D-brane Embeddings 

While the discussion so far was only for the D3/D7 system and the only field strength we 
turned on was Frt, corresponding to finite density in the field theory, the same methods can be 
used to generate analytic embeddings for a variety of probe D-brane systems and worldvolume 
fields. In this subsection we will consider Df^-brane probes in the D3-brane background, with 
q = 5,7. The probe D5-brane wraps AdS^ x S'^ inside AdS^ x S^, and corresponds to adding 
flavor fields confined to a (2+l)-dimensional defect in the SYM theory [19,20]. 

We include F^t and all possible gauge field components Fri with spatial indices i parallel 
to the D3-branes, corresponding to non-vanishing U{1)b current components (Jj) in the field 
theory. We are working at zero temperature, so these currents will not dissipate and can be 
turned on as constant background fields. They do not require an electric field to be turned on. 
For all these cases, the warp factors continue to cancel in the DBI action, and we can obtain 
analytic solutions for the embeddings that basically correspond to Dg-brane embeddings in 
flat space. For the interesting case of non-vanishing electric or magnetic fields, the warp 
factors do not cancel and one has to solve the full-fledged non-linear problem in curved space. 

Our general Dg-brane will be extended along AdSD+2 x 5", so that i runs from one to 
D. Generalizing eq. ( p. 41) , the action for the Dg-brane is 




(2.28) 
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where T^q is the tension of the Dg-brane. Introducing y(r), At{r) and Fri = A[{r)^ the 
Lagrangian for the Dg-brane is 



L = -AAgr"^l + y'2 - (27ra')2Af + {2^a'f ^ Af (2.29) 

where Nq = NfTuqVn with Vn the volume of a unit n-sphere. We will divide both sides of eq. 
( |2.28| ) by the volume of M^'^ so that S^q = / drL. To simplify the formulas let us introduce 
the following notation: 

Y = {„, (2™')A., (2™')A.) , C={c j^y (2.30) 

that is, we construct a vector of fields and the associated conserved charges, with com- 
ponents Ya and Ca- The inner product of these vectors has a Lorentzian signature, i.e. 
diag{+l, —1, +1, +1, +1) for the D7-brane. With this, we have 



L = -^fqr'Wl + Y'^ (2.31) 

and from this we get 

\f r"'Y' 

= ^ ^'^'^ ^" (2.32) 
We algebraically solve for the fields, and find the regulated on-shell action. 



.r2..: r.2 ^ Snq = -Nq (2.33) 



Once again, the lower endpoint of integration depends on the sign of C^. The divergence in 
Soq is cancelled by the action for a single straight Dg-brane with zero field strength. 

So = -Mq r drr^ = -^A"+i, Sren = hm [Soq - So) . (2.34) 
Jo n + 1 A^oo 

{2va')'^ — < for the D7-brane becomes > for the Dg-brane. The Dq^-brane has 

a turn-around point at tq = ^" • The integrals can be done using the machinery in the 



9 



Appendix by changing coordinates to t = rg /r , 

1 /^2N^ ^ f n + l 1 



S„„ = --JV,»{C=)-B -_,-) (2.36) 
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Notice Sren < 0, SO the free energy is positive and hence these embeddings are thermody- 
namically disfavored relative to the embedding with a straight brane/anti-brane pair (which 



still has = 0), just as in section 2.2 



(27rL')^ — > for the D7-brane becomes < for the Dg-brane. The Dg-brane 
does not have a turn- around point and the integrations extend to r = 0. Here the change of 
variables is to u = j^^r'^"', with the result 

1 ,-1, ( My 1 n-l\ , , 

Y. = C.^^N, " \C^\-^B [j^^^d^y ^, ^ j (2.37) 

As in the (^2^1^')'^ — > case for the D7-brane, Sren > 0, so the free energy is negative and 
these embeddings are thermodynamically favored. 



3. The Limit of Zero Mass 

In this section we study the SYM theory at finite temperature and density, and take m = 0. 
At finite temperature, the AdS^ metric in eq. (^) becomes the AdS-Schwarzschild metric, 

fJr^ 

ds'^ = — ^ - f(rQ)dt'^ + r^dx^ + d9^ + sin^ Odsl^ + cos^ Odsls (3.1) 
fire) 

where /(re) = t\ — The black hole horizon is related to the temperature by = ttT. 
The D7-brane embedding is described by 9{rQ). Introducing the worldvolume gauge field 
Frgt = drgAt{rQ), the D7-brane action eq. ( |2.4| ) becomes 



5^7 = -M j dn^rl cos^ eyjl + f{rQ)e'^ -{2^a'YAf (3.2) 

where we have divided by the volume of M^'^, J\f is the same as defined above, and primes 
denote At this point we can plainly see that the system has only one conserved charge, 

d^MrlcQg'e , (^^"O'-^'t (3.3) 

^l + f{r^)e'^-{27:a'YAf 

The algebraic solution for ^^(re) is 



[lira 



(3.4) 



Inserting this into Sdj, we find the on-shell action. 
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/.A 

5^7 = / dr^rl cos^ 9 



^±1M>^ ,3.) 



N'^rl cos6 9 + 



(27ra')^ 



We may obtain the equation of motion for ^(rg) in two ways. We may either derive it from 
eq. (^]^) and then plug in the gauge field eq. ( p.4[ ), or we may eliminate the gauge field at 
the level of the action via a Legendre transform and then derive ^(rej's equation of motion. 
The Legendre-transformed action, Sdi, is 



Sd7 = Sd7- j dr^FrJ-^^ =- j dr^^JTVW^s^M^rl cos^ 9 + (3.6) 

Once again, to cancel the A — > oo divergence in 5^)7 we subtract the action of a single straight 
D7-brane with no gauge field. At finite temperature, the lower endpoint of integration is now 
th-, hence 

So = -M r dr^rl = -jMA^ + jMrjj, (3.7) 

and Sren is defined as above. 

A valid solution of 0(rg)'s equation of motion is simply ^(rg) = 0, corresponding to 
m = and {Om) = in the SYM theory. For this solution, the integrals for At{r) and Sd7 
are identical to those for a black hole embedding with c = 0, eqs. (|]|) and (|l0|), but with 
r — > rg and with the lower endpoint of integration r^- We thus find 




B 



M'^r% 1 1 



(3.8) 

Notice that At{rH) = 0. The value of is given by taking rg oo, for which the first term 
in parentheses goes to (g, For the on-shell action we find 

d2 



4^6 \2TTa' J \ 



Sren = " ^A^-^^ [ -— \ B\ . f ^""^ ; - I • (3.9) 




The T = limit of course reproduces the results of section 2.3 for black hole embeddings with 
c = 0. A straightforward analysis reveals that = —Sren has no meaningful features: in the 
m = limit, the system is stable and exhibits no phase transitions. 

4. Conclusion 

For Nf massive N = 2 supersymmetric hypermultiplet fields at finite U{\)b density coupled 
to the M = A SYM theory in the limits of large- A'^c and large A, we have provided analytic 
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expressions for the hypermultiplets' leading-order (in Nf/Nc) contribution to the free energy 
in two hmits: zero temperature and zero mass. We have found a second order transition 
at zero temperature when fj, = m. In the supergravity description, we have confirmed the 
proposal of refs. [7,18] that, in the canonical ensemble, only black hole D7-brane embeddings 
are physical, while in the grand canonical ensemble black hole embeddings are appropriate 
for fi > m while Minkowski embeddings with constant gauge field (and zero field strength) 
are appropriate when ^ < m. Perhaps the most obvious task for the future is a thorough 
analysis of the phase diagram of this theory in the grand canonical ensemble. 
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Appendix: Beta functions and Incomplete Beta functions 

The integrals needed in the bulk of the text all take the form of Beta functions or incomplete 
Beta functions. For the reader's convinience we collect some of the important properties of 
these functions in this appendix. The Beta function is defined as 



B{a, h) = = f dt{l - tf-H^'^ = r du{l + ur'^^+'^u--^ (1) 

r(a + h) Jo Jo 



and the incomplete Beta function as 

B{x;a,b)= r dt{l - tf-H""-^ = r'^^ ''^(iu(l+n)-('^+*)n''-\ (2) 
JO Jo 

These satisfy the recursion relation 

B{x;a,h) = B{a,h) - B{1- x]b,a) (3) 
and are related to the hyper geometric function as 

B{x; a, b) = a'^xTia, 1 - 6; a + 1; x). (4) 
The divergence in Sdi appears in the expansion, for a = —2/3 about x = 0, 

=-^x-2/3 + 0(xV3) (5) 
where for Soj we will have x ~ A~^ (see for example eq. (pi) below). 



The integrals we need can easily be brought into the form of eq. (^). For notational 
simplicity, we will replace (^2Tra')'^ with just d^. We begin with the — < case of section 
The integral for y{r) and At{r), eq. ( |2.9| ), in this case is 
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with tq = N ^/'^ {(? — d'^)^^^. We change variables to t = r^/r^, 



/<(r) 



6 Jt 



The same change of variables works for the regulated on-shell action eq. ( 2.10| ) 



--Mr 



rg/A6 



dt{l - t)-l/2t-5/3 



-\^rt ( B 



2 1 
3' 2 



2 1 

3' 2 



Now we turn to the d — c > case of section 2.3. The integral for y{r) and At{r 



We change variables to u 



I>{r) 



■2^6 I j2 „2\-1/2 



/>(r) = iAA-V3(d2_c2)-i/3 rdu{l + u)-y\-'/^ 
6 ./n 



2^6 



1 1 



1/3(^2 ^2^ ^ ^l^2^6 + rf2_c2'6'3 



The regulated on-shell action is, with A' = 



'-'1)7 



6 
6 

--M-'/\d^ 



^2)2/3 



dn(l + ^z)-i/V/6 
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